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The  nonlinear  interaction  between  a steady  transonic  mean  flow  and  that 
induced  by  superposed  small -amplitude  oscillations  is  studied  using  the  method 
of  harmonic  expansions.  A truncated  system  of  spectral  equations  is  derived 
which  accounts  for  both  frequency  and  amplitude  dependent  changes  to  the  mean 
flow,  as  would  be  induced  by  Reynolds  stress  feedback,  and  in  particular,  that 
owing  to  the  primary  harmonics.  A three-level  mixed-differencing  scheme, 

conservative  to  leading  order,  is  devised  to  solve  for  the  nonlinearly 

^ 

coupled  mean  and  unsteady  disturbance  flows.  ^Calculations  are  performed 
for  the  NACA  64A006  airfoil  with  an  aft -mounted  quarter  chord  oscillating 
flap,  for  a range  of  subcritical  and  supercritical  Mach  numbers,  reduced 
frequencies,  and  flap  deflection  angles.  The  results  are  compared  with 
existing  linear  and  nonlinear  computations,  and  with  some  experimental  results 
of  Tljdeman;  also,  for  flap-type  oscillations,  the  nonlinear  response  of  the 
mean  flow  to  growing  unsteady  perturbations  can  be  described  in  qualitative 
terms,  and  given  simple  physical  explanation.  At  least  for  the  cases  considered. 
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the  results  of  numerical  experiments  suggest  that  the  nonlinear 
•; back- interaction  on  the  mean  flow,  anticipated  especially  for  low  reduced 
-frequencies,  is  insignificant  even  for  "large  amplitude"  oscillations;  also, 
the  response  of  the  unsteady  loading  to  changes  in  flap  deflection  angle 
^appears  to  be  a linear  one.  The  effects  of  nonlinearity,  however,  are  important 
to  the  extent  that  unrealistic  overpredictions  in  leading  edge  perturbation 
.pressure,  as  obtained  from  linear  finite  difference  calculations,  and  which 
persist  even  at  very  high  subsonic  freestream  Mach  numbers,  are  significantly 
reduced.  The  nonlinear  coupling  therefore  restrains  the  influence  of 
.Jipstream  propagating  or  “receding"  waves  that  do  not  actually  belong  in  the 
..physical  problem.  Calculations  were  carried  up  to  the  point  of  Ehlers' 
..so-called  "Helmholtz  divergence",  beyond  which  uncontrolled  numerical 
.instabilities  were  encountered. -^The  harmonic  method  presented  here  is  not  as 
-^computationally  efficient  as  existing  ADI  methods,  for  example,  the  recent 
-scheme  of  Ball haus  and  Goorjian,  but  it  do^s  yield  some  insight  into  the 
^governing  dynamical  processes,  and  it  may  be  useful  in  nonlinear  flutter 

. analyses.  ^ * ; 

* * • * . .•  ■ / 

1.  INTRODUCTION 

Consider  a two-dimensional  airfoil  in  a subsonic  freestream  executing 
-small -amplitude  harmonic  oscillations  about  a mean  disturbance  flow.  Let 
X»y  and  t be  streamwise,  transverse  and  time  coordinates,  respectively, 
rJ[  be  the  ratio  of  specific  heats,  and  characterize  the  nonoscillating  part 
the  unsteady  configuration  by  a maximum  thickness  h and  a chord  c. 
-Furthermore,  let  , U*,  , ^ and  P denote  the  freestream  values  of 
Mach  number,  streamwise  velocity,  fluid  density,  and  static  pressure.  The 
Steady  transonic  limit  refers  to  TT-^0  and  1,  where  V = h/c  is  the. 

thickness  ratio,  and  the  similarity  parameter  Ks  (1  -M \o)/ is  held 
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fixed.  This  idea  is  extendable  to  unsteady  flows.  For  sufficiently  slow 


temporal  variations,  the  nondimensional  variables  x*  = x/c,  y = T^y/c  and 


t -V  t/c  are  of  order  unity,  and  the  asymptotic  expansion 


$ = +CZ-  (x,y,t)  +...}  for  the  total  velocity  potential  $ yields, 

to  leading  order,  the  perturbation  potential  equation 


{K-a+6pg}?si  -2ffg  - o 


(i) 


In  the  same  low-frequency  approximation,  when  y = TcFfx’.t)  represents  the 
unsteady  airfoil  displacement,  we  have  the  tangency  condition 


Vy(X,°A)’‘  F-  a,i) 


(2) 


■consistently  applied  on  the  slit^  = 0.  The  pressure  coefficient  C 


measures  deviations  of  the  local  static  pressure  P from  , and  is  obtained 

from 

P-7L 
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(3) 


The  Kutta  condition,  applied  on  the  downstream  slit  emanating  from  the 
trailing  edge,  insists  on  pressure  continuity,  and  requires  that  the  jump 


(4) 


Identically  vanish  across  the  wake;  thus,  tv 1 is  a function  of  time 
only,  and  is  instantaneously  equal  to  the  potential  jump  at  the  trailing 
edge.  To  Eqs.  (2)  and  (4)  we  add,  of  course,  regularity  conditions  for  the 


mean  flow,  and  radiation  conditions  for  the  unsteady  disturbance  flow,  but 
In  a form  suitable  for  convenient  numerical  implementation. 

With  airfoil  motions  that  are  harmonic  in  time  and  which,  primarily, 
are  small  unsteady  oscillations  about  a steady  disturbance  flow,  it  is 
plausible  to  linearize  the  unsteady  problem  about  a steady  flow,  where  the 

latter  may  be  taken  to  be  the  solution  for  the  nonoscillating  airfoil.  This 

12 

idea  was  first  given  in  Landahl  ’ although,  at  the  time,  the  extent  to 

which  Landahl's  suggestion  held  true  could  not  be  assessed.  Successful 

.computational  methods  for  mixed-type  boundary  value  problems  admitting 

-Shocked  solutions  have  only  become  available  recently.  In  1971,  Murman  and 
3 

■Cole  successfully  calculated  mixed  flows  with  embedded  shockwaves  using  a 
•type-dependent"  relaxation  method;  their  scheme,  for  nonlifting  flows,  was 

4 

.extended  to  lifting  configurations  in  Krupp  and  Murman  , and  was  further 
Improved  to  allow  for  the  proper  conservative  capture  of  flow  discontinuities 

c 

In  Murman  . These  investigations  have  established  the  computational 
- feasibility  for  practical  transonic  flow  calculations.  The  calculation  of 
'Unsteady  flows,  as  suggested  in  Landahl's  linearized  formulation,  is  also 
-amenable  to  type-differencing  techniques.  Ehlers^  and  Traci,  Albano,  Farr 
and  Cheng^,  for  example,  have  independently  computed  the  unsteady  field 
^-Induced  by  small,  unsteady  oscillations  using,  as  input  to  Landahl's 
.•variable  coefficients,  the  steady  background  flow  as  obtained  from  Murman's 
algorithm  (in  both  of  these  studies,  the  problem  is  Fourier-transformed  in 
^lime,  as  will  be  the  case  in  the  present  paper).  The  foregoing  cited 
^references  form  the  basis  of  modern  transonic  flow  calculations. 

In  linearized  theory,  the  nonlinear  mean  flow  is  obtained  once  from  a 
Murman-Cole  method,  and  its  solution  fixes  the  variable  coefficients  of  a 
linear,  frequency-dependent,  mixed-type  problem.  This  linearization. 


.however,  is  not  strictly  correct:  the  nonharmonic  part  of  the  total 
disturbance  flow  must  change  for  different  disturbance  frequencies  and 
.^amplitudes  because  the  physical  problem  is  inherently  nonlinear.  The  mean 
.flow  cannot  be  regarded  as  completely  fixed  for  all  time;  the  " £ -expansion" 
(that  is,  the  straightforward  linearization)  breaks  down  because  the  nonlinear 
-feedback,  which  arises  from  the  usual  harmonic  interplay,  and  which  builds 
-4jp  over  the  initial  transient  period  during  which  the  oscillations  commence, 

.is  not  described. 

-This  back- interact! on  effect,  generally  assumed  to  be  0(fcz),  where  £. 
-^characterizes  the  smallness  of  the  unsteady  oscillation  amplitude  relative 
;-to,  say,  the  thickness  of  the  mean  geometry,  is  ignored  in  linear  theory. 

*sFor  many  practical  applications,  linear  theory  more  or  less  suffices;  however, 
*?ifor  low-frequency  transonic  flows,  generally  associated  with  large 
-^aerodynamic  loadings  as  well  as  large  unsteady  shock  excursions,  the  neglect 
■=©f  nonlinear  feedback  may  well  be  incorrect.  Thus,  one  might  anticipate 
-iAhat  those  formally  0(S*)  terms  may,  in  fact,  be  large:  they  may  have  a 
significant  effect  ofi  the  actual  mean  shock  location  and  strength,  and 
whence,  on  the  unsteady  perturbation  field  (for  high  frequencies,  the  effect 
^is  probably  less  pronounced,  and  a linearized  approach  should  do).  To 
’"describe  this  back- interact! on,  the  nonlinear  coupling  between  mean  and 
^disturbance  flows  must  be  retained,  and  new  expansions  must  be  introduced 
^=to. reflect  the  parametric  dependence  of  the  actual  flutter  problem  on  • 

-^sthe  freestream  Mach  number,  the  oscillation  frequency  Cd  , and  the  oscillation 
^amplitude  £,  given  some  airfoil  section.  In  essence:  the  £ -expansion 
•vanay  not  retain  its  uniform  validity  over  the  entire  range  of  Cd's.  The 
present  study  explores  the  consequences  of  nonlinearity  for  the  well 
-^documented  NACA  64A006  airfoil.  • 


II.  ANALYSIS 


To  Implement  the  ideas  discussed  in  the  Introduction,  expand  the 
disturbance  potnetial  <p  in  the  real  series 

~ wt  -iJt 

where  1S  the  comPle*  conjugate  of  py  and  higher  harmonic 

terms  are  not  shown.  Linearized  unsteady  theory  defines  a sequence  of 
problems  by  equating  coefficients  of  like  powers  of  £ (the  <pQ  problem,  of 
course,  is  nonlinear).  This,  as  previously  stated,  is  not  always  correct 
because  the  nonharmonic  mean  flow  is  rendered  independent  of  £ and  , 
and  this  is  physically  unrealizable.  To  remedy  this  defect,  a sequential  ordering 
Is  instead  introduced  by  equating  coefficients  of  like  harmonics.  The 
resulting  equations,  if  they  can  be  accurately  solved,  involve  no  approximations 
and  describe  the  physical  problem  exactly.  However,  since  they  are  infinite 
In  number,  some  type  of  truncation  is  needed.  This  is  possible  if  £ ^1; 
for  then,  it  is  likely  that  higher  harmonics,  associated  with  high-order 
amplitude  effects,  can  be  neglected.  We  therefore  retain  in  the  mean  flow 
equation  only  those  "feedback"  terms  arising  from  the  primary  harmonics. 

Hence,  the  expansions  .... 

<v  ^ . iu/'t  . 

= ?.(*,/)  J («) 

~ * £ c°s  $•*>  wt ) 

and 

Ftfjt)  = -ftf)  e } (7) 

e "f( X)  * l ( jrcoS  wt  - <j* si*  u/'t  ) • . 


where  (f>  j = +1  <p ^ and  g = gr  + ig1  , In  the  usual  notation,  lead 

to  three  coupled  nonlinear  boundary  value  problems,  namely,  for  the  mean  flow, 
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for  the  flow, 
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and,  for  the  flow, 

•{K'(ir+0^}  <Py__  -U+0<%.~ +vL  ‘2“%-  < 

X XX  *A  * // 

= ^ on -the  chord  ( 
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The  ideas  underlying  the  foregoing  expansions  were  first  described  in 
Chin®,  where  the  expanded  forms  for  the  corresponding  jump  conditions  may 
be  found.  These  are  suited  for  shock-fitting  analyses,  but  the  overall 
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complexity  of  the  total  formulation  forbids  rigorous  numerical  analysis; 
thus,  a plausible  "shock-capturing"  algorithm,  which  ignores  the  internal 
boundary  defined  by  the  evolving  shockwave,  is  used. 


III.  NUMERICAL  PROCEDURE 

A column  relaxation  scheme  employing  mixed-differencing  and  tridiagonal 
matrices  is  used  which  sweeps  the  computation  box  from  upstream  to 

/»  t 

downstream.  The  field  variables  , p?  and  are  first  initialized  to 
zero.  Eqs.  (8-10)  for  the  mean  potential  $ are  first  solved  by  "sweeping" 
the  flowfield  once,  with  ^ "frozen";  then,  Eqs.  (11-13)  are  solved,  again 

t i 

sweeping  through  once,  with  (p0  and  ^ frozen,  and  finally,  ^ is  solved 
using  Eqs.  (14-16),  with  (po  and  (p^  fixed.  This  triple-sweep  operational 
sequence  defines  one  "iteration".  The  latest  available  values  of  (f0  and 
are  always  used  to  evaluate  the  coefficient  and  "right  side"  arrays, 
effectively  linearizing  the  solution  along  each  line.  It  is  not  clear  how 

to  construct  a strictly  conservative  scheme.  One  needs  to  generalize 

5 • 

Murman's  derivation* for  the  steady  shock-point  operator  for  unsteadiness, 

g 

using  the  somewhat  cumbersome  jump  conditions  in  Chin  , and  this  has  so  far 
proved  unwieldy. 

Conservative  schemes  for  the  steady,  lifting  problem  are  readily 

available.  In  computing  the  mean  flow,  the  0(£?)  right  side  of  Eq.  (8), 

which  represents  unsteady  Reynolds  stress  effects,  is  assumed  to  be  a known 

prescribed  forcing  function.  Murman's  conservative  differencing  is 

q 

Implemented  through  Jameson's  convenient  "switching  function"  formulation. 
For  example,  if  4x  and  Ay  are  constant  nonvariable  mesh  widths  in  the 
streamwise  and  transverse  directions,  respectively,  whose  relative  positions 


r 


are  Indicated  by  the  indices  i and  j,  and  if 
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a fully  conservative  approximation'  to  Eq.  (8),  conserving  mass  across  the 

evolving  shock  is,  with  errors  of  0(£x),  (1  -S.  .)P.  . +S.  , -P.  , . +q.  . 

i,j  i»J  i-l,j  l-l, j Mi,j 

* ^i»j  * w^ere  ^HSi,j  denotes  the  right  hand  side  of  Eq.  (8).  This 
difference  equation  is  solved  through  column  relaxation,  as  suggested  by 
the  regrouping 
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The  airfoil  tangency  condition  and  the  wake  Kutta  condition  are  handled  in 

the  manner  of  Krupp  and  Murman^,  taking  due  account  of  the  RHS.  . array,  and 

■ »J 

the  farfield  boundaries  defined  by  the  computation  box  are  updated  at  the 
end  of  each  sweep.  The  tridiagonal  equations  obtained  along  each  column 
are  linearized  about  "old"  values,  and  the  columns  are  solved  successively 
from  upstream  to  downstream.  The  details  of  this  algorithm  are  straightfowardly 
Imp! emented. 

r t 

For  the  perturbation  flows  ^ and  we  proceed  in  a similar  manner. 

The  construction  of  proper  shock-point  and  parabolic  operators,  in  the  sense 

5 3 

of  Murman  , proved  unwieldy.  Thus,  Murman  and  Cole's  original  type-switching 

was  employed,  and  this  was  implemented  together  with  Ehlers'®  treatment  of 

unsteady  tangency  conditions  and  some  simplified  farfield  boundary  conditions 

of  Seebass  (to  be  described  later).  The  unsteady  solution  for  , therefore. 

Is  not  fully  conservative,  although  the  mean  flow  solution  is,  to  leading 

order.  The  sequential  solution  for  (jP ^ and  then  j completes  the  corrector 

level  of  the  triple-sweep  iteration,  at  which  point  the  mean  potential  is  once 

again  solved.  Both«unsteady  Kutta  conditions  are  implemented  as  in  Ref.  4. 

On  convergence  of  the  relaxation,  calculation  of  the  surface  pressure 

coefficient 

• . . . 

Cf  S - </>l_  Sm  wt  ) (17) 


proceeds  by  streamwise  differentiation  of  linearly  extrapolated  surface 
potentials.  The  definition 
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Is  especially  convenient.  Our  subsequent  results  will  thus  be  expressed  in 

terms  of  the  functions  C = -2  £ , C**  = -2zl3<fir  , and 

P0  r<,x  P] 

r1  B -2  ~CX^(£ 1 . Note  that  £ is  assumed  to  be  proportional  to  the 

\ ** ' 

amplitude  of  the  unsteady  oscillation,  e.g.,  the  flap  deflection,  so  that  the 

f*  1 

functions  and  Cp  , as  defined,  should  remain  nearly  invariant  throughout 
a range  of  sufficiently  small  e ’s,  as  would  be  the  case  in  linear  theory. 


IV.  CALCULATED  RESULTS 

All  of  the  results  reported  here  were  obtained  on  nonvariable  meshes 
consisting  of  80  grids  in  the  streamwise  direction,  with  twenty  over  the 
chord,  and  60  grids  in  the  transverse  direction.  Both  A x and  A'y,  and 
supersonic  and  subsonic  relaxation  parameters,  implemented  as  in  Traci  et  al^, 
were  varied  from  case  to  case  to  provide  stability  and  to  ensure  rapid 
convergence  although,  of  course,  in  runs  designed  to  test  the  effect  of 
changing  a single  parameter,  the  same  convergence  parameters  were  assigned. 

For  simplicity,  in  the  unsteady  problem,  farfield  boundary  conditions  calling 

* 10 
for  = 0 downstream  and  ^ = 0 elsewhere,  as  in  Ballhaus  and  Goorjian 

and  Fung,  Yu  and  Seebass^ , were  used.  The  algorithm  devised  here  is 
expensive  and  slowly  convergent;  typically,  1000  iterations  are  necessary 

' M 

to  produce  residuals  of  0(10  ;,  although  results  sufficiently  accurate  for 
engineering  applications  can  be  obtained,  to  within  10%  error,  in  200  or  so 
Iterations.  The  present  algorithm  also  suffers  from  Ehlers'®  so-called 
•Helmholtz  divergence"  (see  also  Reference  7):  nonconvergence  occurs  for 
sufficiently  high  freestream  Mach  numbers  and/or  reduced  frequencies  k 
based  on  semichord.  Of  course,  the  present  work  concerns  itself  with  low 
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frequency  flows,  and  the  latter  constraint  is  not  very  significant;  however, 
the  former  restriction  on  very  supercritical  flows  prevents  us  from  considering 
those  severe  cases  for  which  large  shock  oscillations  are  anticipated. 
Nevertheless,  the  effects  of  nonlinearity,  even  for  high  subsonic  and  mildly 
supercritical  flows,  are  of  engineering  and  mathematical  interest:  how 
siqnificant  is  the  Reynolds  stress  feedback? 


To  answer  this  question,  the  flow  about  the  NACA  64A006  airfoil  with 
an  aft-mounted  quarter  chord  oscillating  flap  was  studied.  In  particular,  we 
considered  the  parameter  ranges  0.80^  M^-0.85  and  0.00  6 k ^ 0.25  (the 

= 0.85,  k = 0.25  case  was  the  maximum  nondivergent  parameter  set  possible 
for  the  grid  sizes  and  relaxation  factors  used),  and  flap  deflection 
angles  $ = 1°  , 2°  and  3°.  This  run  portfolio  is  especially  useful  since 
comparison  with  existing  linear  results  and  experimental  data  is  possible. 

All  of  the  calculations  reported  in  this  paper  assume  an  u'pitched  mean 
geometry.  We  examine  first  several  results  for  a flap  deflection  of  S - 1°, 
and  a freestream  Mach  number  of  0.80.  The  effects  of  increasing  reduced 
frequency,  shown  in  figures  1,  2 and  3,  correspond  to  k = 0.0,  0.064  and 
0.25.  All  of  the  mean  surface  pressure  coefficients  obtained  are  identical 
to  within  0.5%,  and  fall  between  the  coarse  mesh  and  fine  mesh  results  of 
Traci  et  al.  Good  experimental  correlation  is  also  obtained.  The  unsteady 
force  coefficients,  in  the  convention  of  Ehlers^,  Traci^  and  Tijdeman^’^’^, 

- V*  *?  p j 

are  expressed  in  terms  of  cross-chord  jumps  in  = £Cp*  IS  » that  is. 

In  terms  of  Afc  . In  all  three  figures,  the  hingeline  Ag  £ is  well 
predicted  by  the  nonlinear  method.  Leading  edge  results,  however,  are  far 
from  satisfactory.  For  example,  in  calculating  ^ £ , the  linear  finite 

difference  methods  of  Ehlers  and  Traci  both  appear  nearly  singular  while,  on 


,±he  other  hand,  the  magnitudes  are  underpredicted  by  the  present  theory. 

*Jhe  exact  reasons  for  the  observed  discrepancies  are  unclear.  Certainly,  on 
physical  grounds,  the  effects  of  unsteady  pressures  ahead  of  the  hinge  axis 
-Should  be  minimal  for  high  enough  Mach  numbers;  this  casts  some  doubt  on  the 
results  of  linear  theory,  and  in  this  respect,  nonlinear  theory  may  represent 
.a  needed  improvement.  The  effect  of  leading  edge  bluntness,  too,  was 
--.tested  by  modeling  the  infinite  slope  in  separate  runs  using  different  "large" 
values,  but  the  overall  qualitative  results  for  ^ ^ p’1  were  unchanged  from 
;?those  presented  here.  It  should  be  noted,  however,  that  the  present  results 
-correlated  somewhat  better  with  experimental  data  corresponding  to  higher 
--*Hach  numbers  (all  of  the  calculations  reported  here  assume  an  unbounded 
agfreestream,  whereas  Tijdeman's  experiments  were  conducted  in  a porous  wall 
=arind  tunnel).  Figure  4,  for  Mwr=  0.825  and  k = 0.0,  is  to  be  compared  with 
vFigure  1;  Figure  5 repeats  the  calculation  of  Figure  2,  except  that  0.825. 
—Again,  the  same  qualitative  differences  between  linear  and  nonlinear  results 
-are  evident.  Agreement  between  the  linear  finite  element  results  of  Chan  and 
-JJrashears^  and  the,present  results  in  Figure  2 appears  to  be  fortuitous. 

Figure  6 considers  Mo*=  0.85  and  k = 0.24,  again  with  S = 1°.  The  mean 

^pressure  coefficient  agrees  moderately  well  with  experiment.  The  unsteady 

^peak  pressure  coefficients  are  predicted,  and  the  qualitative  shapes  of  the 

loading  curves  agree  somewhat  with  wind  tunnel  data.  The  experimental  results 

*shown,  as  well  as  the  (dashed)  results  of  an  unsteady  linear  calculation,  are 
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retaken  from  Tijdeman  . 

^Figures  7 and  8 deal  with  the  question  of  increased  flap  deflections. 

£Jhe  results  for  = 0.85,  k = 0.24  and  S = 1°,  given  in  Figure  6,  are 
--Teplotted  on  the  upper  level  of  Figure  7,  here  in  terms  of  C , C*  and 
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Cp  . Below  each  of  these  curves  are  the  corresponding  computed  changes  that 

would  arise  from  increases  in  flap  deflection,  relative  to  the  results  for 

$-  1°.  For  example,  below  the  C curve  for  S = 1°,  the  plots  for 

p0 

, C_  (£=2°)  -C  (S=l°)  and  Cn  (£=3°)  -C  (S=l°)  are  given,  and  so  on. 
p0  p0  p0  p0 

The  curves  for  mean  flow  are  particularly  interesting.  For  flap- type 

oscillations,  an  increase  in  flap  deflection  angle  is  seen  to  induce  an 

underexpansion  upstream  of  midchord  and  downstream  of  the  hingeline,  while 

in  0.5  ^x  ^ 0.75,  approximately,  the  flow  overexpands.  This  same  character 

Is  evident  from  Figure  8,  which  considers  M00=  0.80  and  k = 0.064  (both 

the  reduced  frequency  and  the  freestream  Mach  number  here  are  reduced  relative 

to  their  levels  in  Figure  7).  However,  in  all  of  these  cases,  the  changes  to 

the  C that  would  have  been  obtained  on  the  static  nonosci Hating  airfoil 
P0 

are  probably  insignificant  in  comparison  with  changes,  say,  that  would  arise 
out  of  changes  to  the  computational  parameters.  For  the  unsteady  flow,  the 
3°  corrections  (relative  to  the  £ = 1°  results)  are  approximately  twice 
those  for  £=  2°,  su^esting  near-linearity.  This  is  confirmed  by  the  overall 

smallness  of  the  correction  magnitudes:  changes  in  C„  with  respect  to  changes 

r 1 

In  d are  almost  negligible,  although  in  some  isolated  cases,  5%  differences 

can  be  found  near  the  mean  shock  location  and  near  the  hinge  line.  Note  that 

the  mean  C*  level,  shown  in  Figures  7 and  8,  is  defined  from  K -(K+l)  = 0. 

p0 

Further  increases  in  M,^  or  k beyond  those  values  used  in  Figure  7 gave  rise 
to  uncontrollable  instabilities  of  the  "Helmholtz"  type.  The  question  of 
linearity  will  be  further  reviewed  in  the  next  section,  in  light  of  experimental 
results. 


V.  DISCUSSION  AND  SUMMARY 

When  the  present  study  first  suggested  itself,  the  aim  was  to  demonstrate 

the  predominance  of  nonlinear  feedback,  especially  for  those  critical 

Instances  characterized  by  high  freestream  Mach  numbers,  low  reduced  frequencies, 

or  large  flap  deflection  angles.  At  least  for  the  examples  considered,  however, 

the  back-interaction  on  the  mean  flow  has  been  a weak  effect.  The  foregoing 

results  suggest,  in  fact,  that  linearity  holds  and  remains  applicable  over 

-a  fairly  large  range  of  , k and  £ . This  agrees  with  some  tentative 

■conclusions  of  Fung,  Yu  and  Seebass^.  The  near-linearity  of  the  unsteady 
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.problem  is  also  consistent  with  results  obtained  by  Magnus  and  Yoshihara 
for  the  thicker  NACA  64A410  airfoil,  oscillating  in  pitch  at  Mach  0.72  and 
at  various  frequencies.  These  latter  authors  found  that,  even  for  this 
very  supercritical  flow,  pressures  at  given  locations  and  integrated  quantities, 
such  as  normal  force  and  pitching  moments,  responded  sinusoidally  in  general; 
noticeably  nonsinusoidal  behaviour  in  the  aft  shock  surface  pressure  was 
observed,  however,  for  one  very  low  frequency  case. 
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A second  indication  of  near-linearity  is  seen  from  Tijdeman's 
-experiments.  For  the  NACA  64A006  airfoil  at  zero  pitch  and  zero  mean  flap 
-deflection  angle,  Tijdeman  reports  that  "the  mean  steady  pressure  distributions 
correspond  reasonably  well  with  the  steady  pressure  distributions  ...  obtained 
^on  the  nonoscillating  model"  (these  observations  were  based  on  a flap 
deflection  of  1°,  with  k = 0.25  and  0.80  ^ 0.94).  The  question  of 

linearity,  pursued  in  this  paper,  was  also  explored  by  Tijdeman.  Tests  were 
conducted  to  establish  the  range  of  possible  maximum  flap  deflection  amplitudes 
-wer  which  a linear  relationship  existed  between  S and  the  resulting  unsteady 
loading.  Linearity  was^  established,  but  unfortunately,  due  to  limitations  in 


the  particular  test  set-up,  only  amplitudes  up  to  1°  were  considered. 

On  an  experimental  basis  alone,  Tijdeman's  results  are,  of  course, 
inconclusive:  for  small  deflection  amplitudes,  the  buffer  furnished  by  the 
unsteady  boundary  layer  reduces  flap  effectivity,  so  that  the  extent  of 
"inviscid  nonlinearity"  cannot  be  truly  assessed.  The  calculations  reported 
in  Figures  7 and  8,  however,  do  support  Tijdeman's  conclusion;  these  results, 
purely  inviscid  in  nature,  were  carried  up  to  3°,  and  establish  beyond  doubt 
the  linearity  of  the  unsteady  response.  In  this  light,  linearity  appears 
to  be  meaningful  in  a transonic  field;  the  experiments  of  Tijdeman  and 
Zwaan14,  one  might  further  note,  also  give  a shock  movement  nearly 
proportional  to  the  flap  deflection  for  a sizeable  range  of  flap  angles. 
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This  is  not  to  say  that  nonlinearity  is  unimportant.  Tijdeman  notes 
that  "above  0.875,  the  measured  unsteady  pressures  in  front  of  the 

shock  wave  decrease  strongly,  and  at  Mt->=  0.94,  they  are  zero  ahead  of  the 
hinge  axis."  On  the  other  hand,  both  of  the  existing  linear  finite  difference 
methods®’^  give  contrasting  results.  For  example,  Traci  et  al^  observed 
that  in  both  subcritical  and  supercritical  cases,  the  pressure  perturbations 
over  the  forward  portion  of  the  airfoil  (x  4 0.75)  are  consistently 
overpredicted;  this  is  true  even  in  their  Mo0=  0.94  test  case,  the  calculated 
results  showing  a small  but  almost  uniform  pressure  perturbation  over  the 
forward  airfoil  surface,  whereas  Tijdeman's  results,  which  are  consistent 
with  the  physics,  indicate  negligible  effects.  Traci  attributes  the  observed 
discrepancies  to  real  fluid  effects  and,  perhaps,  to  numerical  inconsistencies 
in  the  implementation  of  farfield  boundary  conditions.  Although  these  are 
distinct  possibilities,  the  global  influence  of  the  nonlinear  coupling, 
retained  in  the  present  analysis,  cannot  be  discounted.  Various  checks  on 
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the  numerical  procedure  have  ruled  out  spurious  effects  tied  to,  say, 
coding  technique  or  leading  edge  modeling.  Thus,  it  seems  that  the  added 
nonlinear  coupling  serves  to  restrain  the  influences  of  upstream-propagating 
or  "receding"  waves  when  they  truly  do  not  belong  in  the  physical  picture. 

This  restraining  effect,  as  is  apparent  from  Figures  1-5,  reduces  the 
perturbation  pressure  level  of  linear  theory,  and  leads  to  somewhat  of  an 
underprediction.  Much  of  this  remains  speculative;  for  example,  the 
reasons  for  the  differences  between  linear  finite  difference  and  linear  finite 
element  models  are  still  unclear.  But  it  may  be  that,  although  the  unsteady 
response  to  changes  in  deflection  amplitude  is  linear,  the  exact  pressure 
magnitudes  must  be  obtained  using  a nonlinear  theory.  A similar  situation 
arises  in  Whitham's  sonic  boom  theory;  there,  the  linear  equation  is  valid, 
provided  the  "nonlinearly  corrected"  characteristics  are  used  in  the  farfield. 

In  obtaining  the  results  of  Figures  7 and  8,  extreme  caution  was 
exercised  to  assure  that  calculated  differences  in  surface  pressure  for  various 

‘deflection  angles,  which  are  small,  are  larger  than  changes  that  would  be 

% 

obtained,  say,  from  iteration  to  iteration.  For  consistency,  in  each  of 
Figures  7 and  8,  results  are  presented  at  the  same  converged  iteration  number. 
The  near-constancy  of  the  mean  flow  to  changes  in  S and  k,  as  obtained  in 
these  runs,  does  not,  of  course,  imply  that  this  is  always  the  case;  in  fact, 
the  time-marching  ADI  scheme  of  Ballhaus  and  Goorjian™,  using  the  small 
disturbance  formulation,  has  successfully  modeled  all  of  Tijdeman's  unsteady 
flow  configurations,  namely,  the  so-called  Type  A,  B and  C shock  motions.  It 
does  suggest,  however,  that  truly  nonlinear  effects,  at  least  for  the 
NACA  64A006  airfoil,  are  probably  confined  to  high  subsonic  Mach  numbers 
very,  very  close  to  unity.  . . 


No  attempt  has  been  made  at  obtaining  quantitative  agreement  with 

experiment.  Experimental  uncertainties  (e.g.,  shockwave-boundary  layer 

interaction,  local  separation,  turbulence  and  wall  interference),  together 

with  computational  shortcomings,  especially  in  the  numerical  treatment  of 

leading  edge  expansions  over  blunt  noses,  make  any  "agreement"  fortuitous. 

Accurate  and  meaningful  comparisons  with  other  numerical  work  are  also 

difficult;  differences  in  mesh  definition,  choice  of  relaxation  parameters, 

implementation  of  farfield  boundary  conditions,  and  specification  of  mean 

flowfields  (e.g.,  both  References  6 and  7 use  a nonconservative  ^(x'.y*)) 

contribute  to  the  basic  problem.  Perhaps  the  most  pressing  issue  is  the  need  to 

model  real  fluid  effects  accurately.  Some  headway  has  been  made  in  the 
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pioneering  work  of  Magnus  and  Yoshihara  , however,  but  clearly,  much  more 
work  needs  to  be  done. 

Although  the  nonlinear  harmonic  method  appears  to  give  improved  results 
over  a wider  range  of  flow  parameters,  a practical  limitation  on  the  approach 
is  Ehlers1®  "Helmholtz  divergence"  boundary,  beyond  which  uncontrollable 
numerical  instabilities  arise.  Until  this  problem  is  overcome,  harmonic 
approaches  are  likely  to  be  of  limited  usefulness  in  transonic  computations. 

Of  course,  the  present  method  was  never  intended  for  routine  use,  but  rather, 
to  assess  the  extent  to  which  linear  theory  can  model  a non! inear  problem. 
Efficient  time-marching  schemes,  such  as  the  finite  difference  ADI  method  of 
Ballhaus  and  Goorjian^,  treat  the  nonlinear  shock  motion  correctly  and 
account  for  higher  harmonic  generation.  The  present  work  points  to  the 
need  to  develope  accurate  nonlinear  models. 
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FIGURE  2/  MACH  0.80,  REDUCED  FREQUENCY  0.064,  FLAP  DEFLECTION  1°  CASE. 


FIGURE  3.  .MACH  0.80,  REDUCED  FREQUENCY  0.25,  FLAP  DEFLECTION  1°  CASE 
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FIGURE  6.  MACH  0.85,  RE0UCE0  FREQUENCY  0.24, 
FLAP  DEFLECTION  1°  CASE. 


• (NONLINEAR  HARMONIC  METHOD) 


FIGURE  8.  MACH  0.80,  REDUCED- FREQUENCY  0.064,  FLAP  DEFLECTIONS  1°,  2°  AND  3°  CASES. 
(NONLINEAR  HARMONIC  METHOD) 


